Large-scale group-hierarchical DEMATEL method for complex systems

Existing Decision-Making Trial and Evaluation Laboratory (DEMATEL) methods are mostly suitable for simple systems with fewer factors, and lack effective integration of expert knowledge and experience from large-scale group populations, resulting in a potential compromise of the quality of the initial direct relation (IDR) matrix. To make DEMATEL better suited for the identification of critical factors in complex systems, this paper proposes a hierarchical DEMATEL method for large-scale group decision-making. Considering the limitations of expert knowledge and experience, a method based on expert consistency network for constructing the expert weight matrix is designed. The expert consistency network is constructed for different elements, and the weights of experts in different elements are determined using the clustering coefficient. Following the principles of the classic DEMATEL method, the steps for identifying key elements in complex systems using the large-scale group-hierarchical DEMATEL method are summarized. To objectively test the effectiveness and superiority of the decision algorithm, the robustness of the algorithm is analyzed in an interference environment. Finally, the superiority of the proposed method and algorithm is verified through a case study, which demonstrating that the proposed decision-making method is suitable for group decision-making in complex systems, with high algorithm stability and low algorithm deviation.


Introduction
In the knowledge economy era, the number of influence factors for complex organizations and systems is increasing, and the complex interrelationships between factors lead to complexity and uncertainty in management of complex organizations and systems.Accurately identifying key factors and clarifying the importance between factors have become important research topics in organizational management and decision-making, leading to the emergence of a series of evaluation and decision-making methods [1][2][3].
Professors Gabus and Fontela proposed the Decision-Making Trial and Evaluation Laboratory (DEMATEL) algorithm [14] in the 1970s.The algorithm identifies key elements by constructing a complex interrelationship diagram between factors and is used for analyzing complex systems.
Compared to the other methods, the DEMATEL method, based on graph theory, provides a more intuitive representation of the complex interrelationships between factors and its calculation process is more simple and straightforward.It has been widely promoted and applied in many fields such as system engineering, and management science.Costa Federica et al. investigated the role of human factors in promoting the establishment of sustainable continuous improvement (SCI) environment [15] by the DEMATEL, Huang et al. analyzed the key elements of circular supply chain management (CSCM) [16], Shahriar et al. analyzed the complex influence factors in the Covid-19 vaccine supply chain [17].In addition, there are many application cases, which are not listed one by one here.Furthermore, many scholars have improved the traditional DEMATEL method.Some researchers are committed to combining the DEMATEL with other methods.Mohammad et al. improved DEMATEL by combining it with the best-worst method (BWM) and Bayesian network (BN), and applied it to safety management [18].Sirous et al. studied the problem of selecting technical suppliers by constructing the Delphi-DEMATEL-ELECTRE method [19].In addition to considering the combination with other methods, some researchers have improved the DEMATEL scale method and expanded the expression of expert judgment information, mainly about point estimation judgment information [20,21], fuzzy number estimation information [22][23][24][25], and grey number estimation information [26][27][28].In addition, some researchers have studied the normalization processing of the decision matrix in DEMATEL to solve the problem that the normalized matrix does not converge in some cases.Michnik et al. proposed a new DIM normalization processing method to solve the problem of the unsolvable TIM matrix [29], and Chen et al. expanded the normalization method for the influence matrix in DEMATEL [30].
However, both DEMATEL method and its improved version require experts to analyze each factor in pairs.In a complex system with n factors, experts need to make n(n−1) judgments.When n is larger than 10, this will require a significant workload, leading to emotional fatigue and boredom that may affect the decision-making results.This limitation restricted the application of DEMATEL in complex systems with numerous factors, and it was not given enough attention by researchers.In 2021, Du proposed a hierarchical DEMATEL method [31], which adopts hierarchical decomposition to divide complex systems into several subsystems.Experts only need to make judgments on subsystems, and then use the proposed method for data integration.This method greatly simplifies the DEMATEL calculation steps and processes in complex systems, reduces the workload of experts, and has attracted the attention of many scholars.Moreover, considering the limitations of individual expert knowledge and experience, Du further introduced group decision-making based on the hierarchical DEMA-TEL [32], however, the group decision-making method is only suitable for small expert groups.When dealing with complex system identification, a large-scale expert group with not less than 20 people is usually needed [33,34].In addition, when a large-scale expert group makes decisions, the rationality of expert weights directly influences the accuracy of decision-making results.However, the weight of experts in Du's study is subjective, lacking persuasive expert weight calculation.The references mentioned above were subjected to comparative analysis, as shown in Table 1.
In summary, although hierarchical DEMATEL has improved the shortcomings of traditional DEMATEL method and can effectively analyze factors in complex systems, how to integrate and utilize the wisdom of large-scale expert groups and improve the scientificity of the IDR matrix of hierarchical DEMATEL method is a new and urgent problem to be solved.This article believes that as the number of factors in complex systems increases and the relationships between factors become more complex, decision making by a small group of experts may not be sufficient to cope with such complexity.It is necessary to determine the direct influence matrix in the hierarchical DEMATEL method by a large group of experts.Decision making by a large group of experts presents the following characteristics: 1.The group size is relatively large, usually consisting of no fewer than 20 decision-making experts; 2. The decision-making problem exhibits multidimensional, complex, and stochastic attributes; 3. High consistency requirements need to be met among the group.
When solving problems of large-scale group decision-making, the main difficulties are as follows [35]: 1.There are significant differences among decision-makers.It is necessary to identify the status of each decision-maker and assign corresponding weights to achieve scientific evaluation results.
2. Due to the large size of the decision-makers, it is important to use effective methods to gather the opinions of large groups to avoid leverage effects caused by intentional praise or criticism during the evaluation process.
3. When group opinions are relatively scattered, it is necessary to effectively coordinate the differences in preferences among decision-makers to maximize the satisfaction of largescale group decisions.Identifying the impact of risk factors and sources of information on the decision-making process Methods composition application [19,28] Gray DEMATEL Studying the causal relationships of influencing factors in the decision-making process Methods composition application [22,25] Fuzzy DEMATEL Estimate and map the suitability classes of ecotourism potentials in the study area of "Dunayski kljuc" region (Serbia); Analyzing the facilitating factors for supply chain responsiveness Methods composition application [26] Gray DEMATEL +ANP Explores favorable methods to evaluate the green mining performance (GMP) of underground gold mines Methods composition application [30] DEMATEL A new matrix normalization method has been researched and proposed Innovation in Method [31,32] Hierarchical DEMATEL The hierarchical DEMATEL method has been proposed to make the DEMATEL method applicable to complex systems with many factors; based on the proposed hierarchical DEMATEL method, a program for small-group experts to reach consensus has been designed Innovation in Method https://doi.org/10.1371/journal.pone.0288326.t001 In the problem of large-scale group decision-making, how to objectively determine the weight of each expert is a key issue.However, this issue is often overlooked by researchers.Chen pointed out that only 41% of the cited papers in group decision-making problems mentioned the determination of expert weights [36].
There are three main methods for determining expert weights in existing research: subjective method, objective method, and comprehensive method.For the subjective method, expert weights are calculated based on factors such as age, attitude, and experience, and the mutual evaluation of experts [37].Multiplicative Analytic Hierarchy Process (MAHP) [38], Simple Multi-Attribute Rating Technique (SMART) [39], and Delphi [40] are key methods for subjective expert weight determination.The objective method is based on the evaluation performance of experts, using individual decision matrices (IDMs) proposed by each expert as the main basis for judgment, and assigning different weights, which is usually more objective [41].The weight is mainly determined based on the degree of closeness between the expert's individual decision and the group decision [42][43][44].In addition, some researchers have constructed expert opinion adjustment mechanisms, attempting to achieve consistency of group opinions by adjusting experts' weights or decisions as much as possible.Pang developed a nonlinear programming model and determined expert weights by maximizing group consensus [45], where expert weights were adaptively adjusted based on their decisions.Yang used a fixed-point iteration method to adjust expert weights multiple times [46].
However, from the relevant research on expert weights, it can be seen that when the research problem is a multi-attribute decision-making problem, almost all methods use a weight value to represent the expert's evaluation status for all attributes.For example, literatures [46][47][48] all use a weight to represent the expert's performance in all fields.
In reality, each expert has limited knowledge and abilities in their own expertise.Using a single weight to represent an expert's status in all fields is unreasonable, which cannot fully reflect the expert's professionalism to highlight their important position in their research field.Repeatedly iterating to solve expert weights to achieve consensus of group opinions often puts pressure on experts who adhere to their own opinions, forcing them to give up their decisions and ideas.
Finally, since decision-making problems are subjective progresses, the results of decisions do not have a correct answer, only subjective judgments of "reasonable" or "unreasonable".It's hard to judge whether one's method is superior to others.Most studies usually compare the decision results of their proposed method with other methods through numerical calculation examples.When inconsistencies occur, researchers often use subjective analysis to explain the rationality of their method, which often lacks persuasiveness.How to use a verification experiment to demonstrate the superiority of the proposed method instead of subjective analysis and judgment is an aspect that almost no researchers have studied.
Therefore, based on the above analysis, we find that the existing research on using DEMA-TEL method to identify factors in complex systems has the following shortcomings: 1.The traditional DEMATEL method can only be applied to situations with fewer elements.
When the number of elements n>10, it will significantly increase the number of expert judgments and workload.
2. Although Du attempted to solve problem (1) by using group hierarchical DEMATEL method, the expert group was small and could not better reflect the wisdom of the expert group.In addition, the expert weight was subjectively given and lacked convincing objective calculation.
3. When calculating expert weights using the objective method, the importance of experts in all attribute fields is often measured by a single weight value.In reality, experts often have a certain disciplinary background, and the degree of specialization of their judgment may be higher for some factors in complex systems, but lower for other factors.Using a single weight to determine the expert's judgment status for all factors is not scientific, and a more targeted method should be adopted.
4. There is a lack of convincing means to test the effectiveness and superiority of decision-making methods.Almost all studies subjectively analyze the differences in decision results between different methods when demonstrating the superiority of their methods.This inevitably leads authors to analyze in a way that favors their own methods and lacks persuasiveness.
Based on this, the author proposes a hierarchical DEMATEL method for large-scale group decision-making to identify key factors in complex systems and address the following issues: 1.In response to the heavy workload for experts in identifying complex systems using traditional DEMATEL methods, the hierarchical DEMATEL method proposed by Du is used to reduce the workload of experts in identifying elements of complex systems.The hierarchical DEMATEL method is combined with the large-scale group decision-making, and the expert number is not less than 20, which improves the quality of decision-making.
2. In response to the problem of expert weight solving, considering the potential influences of experts' knowledge, background, and profession, the factor setting weights are distinguished to construct the weight matrix of expert decision-making.Based on the IDR matrix of experts, the expert consistency network for a certain factor is constructed based on the performance of experts in scoring this factor, combined with the clustering coefficient of the weighted network to represent the consistency of experts.The weight of experts in scoring this factor is determined by their contribution to the clustering coefficient of the weighted network, and the weight matrix of expert decision-making is formed for all factors.This method avoids measuring the performance of experts in scoring in all fields with a single weight value, and can well represent the consistency of the group.
3. In order to objectively demonstrate the effectiveness of the proposed method, interference scenarios are set to analyze the robustness of the decision-making algorithm, and subjective analysis of the decision-making results is avoided.The stability and deviation of the decision-making method after interference are analyzed.The interference scenario refers to the implementation of interference on the original expert data of a certain scale to simulate the judgment deviation of experts.The stability index refers to the degree of change in the decision-making result after random interference, and the deviation index refers to the degree of deviation of the decision-making result after interference from the true value.Obviously, when the stability of the algorithm is high, it indicates that the decision-making algorithm will not easily change the decision-making result due to the influence of disturbance.When the deviation of the algorithm is low, it indicates that the decision-making algorithm can ensure a result closer to the true decision even if it is interfered with.
The innovations of this article are: 1. Introducing a new method for identifying the weights of experts in large-scale groups.This method assigns different weights to different indicators, abandoning the practice of using a single weight value to represent the decision-making status of experts under all indicators, in order to address the unique characteristics of each expert in terms of knowledge, skills, experience, and personality.
2. Using the network clustering coefficient to describe the consistency of expert groups in scoring the same indicator, and calculating expert weights through the consistency between experts and the group, to maximize the requirement of opinion consistency in large-scale group decision-making.
3. The methods involved are more suitable for analyzing the correlation between various factors within complex systems and identifying key factors.It can not only reduce the workload of experts but also improve the scientificity of decision-making results.
4. Instead of analyzing the effectiveness of decision-making algorithms through subjective methods as in other studies, this article constructs interference scenarios to analyze the stability and bias of algorithm results when expert decision-making data is interfered with, which is more convincing.
The rest of this study is organized as follows.Section Preliminaries introduces the basic knowledge, mainly including the introduction of traditional DEMATEL method and hierarchical DEMATEL method.Section The Proposed method mainly including the construction of expert consistency network, calculation of weighted network clustering coefficient, calculation of expert weight matrix, collection of expert opinions, and overall calculation steps.In Section Case presentation and Methodology analysis, stability and deviation indicators of the algorithm are introduced, and numerical calculations and comparative analysis of different methods are performed.In Section Conclusion, our conclusions, contributions and innovations are explained.

Preliminaries
This section introduces the DEMATLE method and the hierarchical DEMATEL method.The hierarchical DEMATLEL method is a new method based on the DEMATEL method proposed in [31], which decomposes the complex system into several subsystems, invites experts to score the degree of influence between elements within each subsystem, and finally turn the set of IDR matrices of all subsystems into a super IDR matrix.This method can effectively reduce the workload of experts and consider the hierarchical characteristics of complex systems.Each of these two methods is described below.

DEMATEL method
The DEMATEL method is a structural model expansion method used to establish and analyze the interactions between complex criteria and oriented to factor analysis of complex systems, the basic elements of the DEMATEL method are as follows.
1. Determine the IDR matrix between the elements Suppose a system F contains N elements, denoted as F = {f 1 ,f 2 ,� � �f N }, and experts are invited to judge the degree of direct influence among these N elements using a scale of {0,1,2,3,4}, representing "no influence", "low influence", "medium influence", "strong influence", and "very strong influence", respectively.The degree of influence of element f i on element f j is recorded as x ij 2 f0; 1; 2; 3; 4g; i; j ¼ 1; 2; . . .; N, the IDR matrix X = [x ij ] N×N is constructed according to x ij , and when i = j, x ij = 0, it represents no influence of the element itself.

IDR matrix normalization
Normalize the IDR matrix constructed by experts [49] ε is a non-Archimedean infinitesimal, the role of ε is to ensure that the infinite powers of the normalized IDR matrix can converge to zero in order to satisfy the conditions for the subsequent third step of constructing the comprehensive influence matrix.

Constructing the comprehensive influence matrix
The comprehensive influence matrix T is 4. Calculation of causality and centrality of each factor Calculate r i = ∑ j t ij and d j = ∑ i t ij based on the comprehensive influence matrix T. The sum of each row element of matrix T, denoted by r i , represents the sum of the influence of factor f i on other factors, which is called the influence degree of factor f i .The sum of each column, denoted by d j , represents the sum of the influence of other factors on factor f i , which is called the being influenced degree of factor f i .Let r i +d i be the centrality of f i , which characterizes the relative importance of factor f i in the system.Let r i −d i be the causality of factor then f i is a receive factor.

Calculate the weights of each factor
The weight of factor f i is ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi In summary, the traditional DEMATEL method can be summarized in the following steps.
Step 1, experts are invited to make decisions on the system elements and construct the IDR matrix X;.
Step 2, the IDR matrix is normalized, and the normalization matrix H is obtained by Eq (2); Step 3, the comprehensive influence matrix T is constructed by Eq (3); Step 4, the centrality and causality of the elements are calculated from the comprehensive influence matrix T.
Step 5, the relative importance of each factor is calculated by Eq (4).

Hierarchical DEMATEL method
The traditional DEMATEL method requires the expert to compare each element pairwise, which is suitable for simple systems with few elements.However, when there are many elements in the system, determining the IDR matrix requires a huge amount of work (if there are n elements in the system, experts need to make n(n−1) judgments).This can easily lead to expert mental fatigue and boredom.In addition, complex systems generally have hierarchical characteristics, which cannot be reflected in the traditional DEMATEL method.To address these issues, Du proposed the hierarchical DEMATEL method in reference [31], which is suitable for identifying key elements in complex systems that contain many system factors and have hierarchical characteristics among them.The method mainly includes the following contents: 1. Hierarchical decomposition of the system Hierarchical decomposition mainly includes vertical decomposition and horizontal decomposition.Horizontal decomposition focuses on dividing the critical factor identification problem of complex systems into several simple problems and vertical decomposition focuses on dividing the complex system into multi-level subsystems under a specific rule.
Horizontal decomposition provides the rules for making vertical decomposition [31].
The complex system F is decomposed into subsystems according to the horizontal and vertical decomposition.As shown in Fig 1, the complex system F can be decomposed into subsystems F 1 ~FQ , and F 1 ~FQ can be decomposed into subsystems the decomposition stops when the system F is decomposed downward to the factors

The IDR matrix of subsystem
The experts are arranged to score all the subsystems, such as F, F 1 ~FQ and f 1 � f 1 N 1 , to obtain the IDR matrix of each system.For example, when system F contains Q subsystems, the IDR matrix of system F is denoted as X = [x qq 0 ] Q×Q , x qq 0 refers to the degree of direct influence of the q subsystem in F on the q 0 subsystem; similarly, the IDR matrix of system F q containing N q subsystems can be denoted as X q ¼ ½x q nn 0 � N q �N q , and the superscript q of x q nn 0 is used to denote the system to which it belongs.

Calculate the super IDR matrix of the total system
The super IDR matrix of the total system F is obtained by integrating all the subsystem direct influence matrices after correction, which indicates the degree of direct influence among all elements.The integration rules are where the calculation rules for the elements are When q = q 0 , which means that the two subsystems are identical, x qq 0 ij ¼ x q ij , the elements can be obtained directly from the expert scoring matrix for that system; When q6 ¼q 0 , it means that the two subsystems are different, x qq 0 represents the degree of direct influence of the subsystems F q and F q 0 , the subscript number represents the order of the subsystems F q and F q 0 in the IDR matrix of their superior system.z q i represents the centrality of factor f q i in subsystem F q , i.e., ii 0 and t q ii 0 are the elements in the comprehensive influence matrix T q of the subsystem F q with T q ¼ ½t q ii 0 � N q �N q and T q is normalized according to the IDR matrix X q of the subsystem F q using steps (2) to (3).For the convenience of example, the above describes the simple case of two levels.When the level decomposition of a complex system involves multiple levels, the modified IDR matrix of the subsystem needs to be derived sequentially from the bottom to the top level, and the recursive integration from low to high forms the super IDR matrix, and the specific process is detailed in the literature [31].

Calculation of elemental importance and centrality
� X is brought into the DEMATEL method as the IDR matrix, as in steps (1) to (5) in Section DEMATEL method, and the importance and centrality of each element are calculated.

The proposed method
In this section, we introduce the hierarchical DEMATEL method to large-scale group expert decision making.We will construct an expert consistency network based on the performance of the large-scale group experts when scoring the same factors, express the consistency of the experts through the weighted network clustering coefficients, determine the expert weights using the contribution of each expert to the consistency, each expert will have different weights when scoring different factors, express the different abilities shown by the experts in different fields through the weight matrix, and finally the resulting weighted IDR matrix of each system.The weighted IDR matrix will be integrated using the hierarchical DEMATEL method to obtain the centrality and the causality of each factor, and realize the large-scale group hierarchical DEMATEL decision.On the one hand, this paper combines large-scale group experts with the hierarchical DEMATEL, proposes a new method of pooling group wisdom, and improves the authority and scientificity of the IDR matrix of subsystems; on the other hand, it integrates the ability of different experts in different fields, and highlights the professional status of experts in a certain field, because even some authoritative experts, who make judgments in certain fields, are not as scientific as experts who are good at that field.

Problem description and hypothesis
The hierarchical DEMATEL method is used to make decisions on the importance of various elements within a complex system.Firstly, the complex system is divided into several subsystems according to levels and categories, and then m experts are organized to judge the degree of influence among the subsystem factors to derive the IDR matrix.Due to the large number of levels and subsystems, the information of the parent system at each level is used to name a certain subsystem in the form of F q 1 �����q pÀ 1 , and q 1 . ..q p−1 represents the information of the parent system at each level.
For example, a subsystem with subscript q 1 �q 2 , q 2 means it's serial number in its parent system, q 1 is the serial number for its parent system at the higher level, and so on.In Fig 2, system F 1�1 represents the subsystem of the second level, and its parent system is the first system of the first level.The IDR matrix given by the nth expert to the subsystem F q 1 �����q pÀ 1 is denoted as X ðnÞ q 1 �����q pÀ 1 , which includes pairwise comparisons of the degree of influence between the factors involved, using a scale of {0,1,2,3,4}, representing "no influence", "low influence", "medium influence", "strong influence", and "very strong influence".The interrelationships between the systems are shown in the Fig 2.
The problem to be solved in this paper is: How to fully explore the IDR matrix of m experts on the subsystem analysis, determine the weight of each expert in scoring different factors, and obtain the weight matrix of experts to calculate the weighted IDR matrix � X q 1 �����q pÀ 1 of the subsystem, and finally achieve the identification of the importance of each element for all factors based on the hierarchical DEMATEL method.

Building expert consistency networks
Suppose that m experts are organized to score the subsystem F q 1 �����q pÀ 1 , which has K elements, where the IDR matrix given by the nth expert is denoted as X ðnÞ q 1 �����q pÀ 1 2 R K�K , X ðnÞ q 1 �����q pÀ 1 ¼ x The element x ðnÞ ij represents the influence degree of factor i relative to factor j in system F q 1 �����q pÀ 1 made by the nth expert, the process of constructing the IDR matrix by the organization experts is shown in Fig 3.
Due to the differences in professional background and competence knowledge of experts, certain subjectivity and deviation will occur in scoring decisions, and different weights are needed to measure the importance of experts' decisions.In this paper, we use a weight matrix to determine the importance of experts for different factors instead of a single weight value, and establish an expert consistency network through the scoring performance of each expert to calculate the weight of experts for that factor.
Definition 1: Consistency of experts Existing experts a and b make judgments on the degree of influence of factor i on factor j in system F q 1 �����q pÀ 1 .x ðaÞ ij is the decision value of expert a and x ðbÞ ij is the decision value of expert b.Define the degree of agreement between experts a and b as Where a,b2{1,2,� � �m}, T N is the difference between the maximum and minimum values in the DEMATLE evaluation scale, and in scales {0,1,2,3,4} of this paper, T N = 4.
The rationality and nature of the above definition is discussed as follows: Property 1: σ ab 2[0,1], when rating the importance of factor i compared to factor j, if expert a and b give equal ratings, i.e., when x ðaÞ ij ¼ x ðbÞ ij , σ ab = 1, indicating that the two experts' agreement degree to reach the maximum value; when the difference between the ratings is the largest, i.e., when jx ðaÞ ij À x ðbÞ ij j ¼ T N , σ ab = 0, indicating that the two experts' agreement degree reach the minimum value.

Property 2:
As the ratings given by the two experts become closer, the value of σ ab will be bigger, indicating that the agreement degree between two experts is also greater, and σ ab = σ ba , i.e., the relationship between the two decision makers is symmetric.
When scoring element x ij in the system F q 1 �����q pÀ 1 , an undirected weighted network is constructed based on the consistency exhibited by m experts.The m experts are the nodes in the complex network, and the agreement degree among experts is the weight of the edges.Assuming that m experts score element x ij as x ð1Þ ij ; x ð2Þ ij ; � � � x ðmÞ ij , the consistency matrix for element x ij formed by m experts is The matrix Y q 1 �����q pÀ 1 ij can be considered as the consistency network adjacency matrix when m experts score the element x ij in the system F q 1 �����q pÀ 1 .The network can be denoted as represents the set of m experts as nodes, E q 1 �����q pÀ 1 ij ¼ fe 1 ; e 2 ; � � � e h g represents the set of edges, and the weights of the edges are the agreement degree σ ab among the experts.
Similarly, the expert consistency network for all other elements can be constructed based on the expert scoring performance.Obviously, for subsystem F q 1 �����q pÀ 1 , since it contains K elements, K×K consistency networks can be formed, and each network corresponds to an adjacency matrix, and the adjacency matrices of all networks will form the Super Consistency Matrix.The process of constructing an expert consistency network is shown in Fig 4.
The later section will effectively reveal the intrinsic connections among the expert members by further analyzing the consistency network of experts in order to clarify the weight of the experts' scores for each element.

Solving the weighted IDR matrix based on the weighted network clustering coefficient
After constructing experts consistency networks, we can study the relationship between experts based on the network features.The Clustering Coefficient of a network is the ratio of the connectivity between any two nodes in the network to the connectivity between the neighboring nodes they share.
In general, suppose a node v i has k i edges connecting it to other nodes, and these k i nodes are called neighbors of node v i .Obviously, there are at most C 2 k i possible edges between these k i nodes.The ratio of the actual number of edges E i and the total number of possible edges C 2 k i between k i neighboring nodes of node v i is defined as the clustering coefficient C i of node v i , i.e.
In simple terms, it is the ratio of the actual number of connections around a node to the theoretical maximum number of connections.The clustering coefficient C of the whole network is the average of the clustering coefficients C i of all nodes v i .That is Clearly, when all nodes are isolated and there are no connecting edges, C = 0; and C = 1 only when the network is globally coupled, with any two nodes directly connected.
The clustering coefficient represents the degree of closeness and stability of groups formed in the network.When the clustering coefficient is higher, it indicates that the neighbors of the node are closer and the resulting clustering groups are more stable.Since this paper constructs a network based on the scores given by experts for a certain factor, it is only necessary to determine whether the clustering groups for the same factor judgment scenario are stable.The property of the clustering coefficient can just reflect the consistency of the expert group.The expert consistency network in this paper is a weighted network.Onnela [50] studied the clustering coefficient of weighted networks, and defined the clustering coefficient of nodes v i in a weighted network as: In this equation, w ij represents the edge weight between node v i and v j , wij ¼ w ij =maxðw ij Þ is the normalized weight, and k i is the degree of node v i .Combining Eqs (10) and (13), when all experts score the elements x ij in subsystem F q 1 �����q pÀ 1 , the clustering coefficient vector can be obtained based on the consis- Where We can perform exponential normalization on Eq (14) with the exponent m being the number of experts, to obtain the weight value of the nth expert in network G q 1 �����q pÀ 1 ij .Normalize Formula ( 14) by exponentiation, where the exponent m is the number of experts, to obtain the weight value w n of the nth expert in network G q 1 �����q pÀ 1 ij , there is To distinguish the subsystem information and elements targeted by the weight values, we add a superscript q 1 �� � ��q p−1 :ij to the weight values in Eq (17), so that the weight value of the nth expert when judging the degree of influence of factor i on factor j in system F q 1 �����q pÀ 1 can be expressed as w q 1 �����q pÀ 1 :ij n , where Repeating Eqs (11) to (18), we can obtain the weight matrix of experts for all factors in the system F q 1 �����q pÀ 1 .The weight matrix of the nth expert for the system F q 1 �����q pÀ 1 is By combining the weight matrix and the IDR matrix, the weighted IDR matrix � X q 1 �����q pÀ 1 for subsystem F q 1 �����q pÀ 1 is Based on Eqs ( 8) to (20), we can obtain the weighted IDR matrix for each subsystem by constructing an expert consensus network, calculating expert weights based on the weighted network clustering coefficient, and synthesizing the weighted IDR matrices for each subsystem.By using the weighted IDR matrices for each subsystem as the IDR matrix of the hierarchical DEMATEL method, we can analyze the causality and centrality of all elements in a complex system by the hierarchical DEMATEL method.
In summary, the main process of the proposed method in this article is shown in Fig 5.

Calculation steps of the proposed method
Based on the previous description, the computational steps of the large-scale group decision hierarchical DEMATEL method proposed in this paper are sorted out as follows.
Step 1: To classify a complex system hierarchically For a system F = {f i |i = 1,2,. ..,y}, the influence factors f i are divided into subsystems by hierarchy and by attributes, and the subsystems are denoted as F q 1 �����q pÀ 1 ; Step 2: Constructing the IDR matrix for each subsystem.The IDR matrix obtained from the judgment of subsystem F q 1 �����q pÀ 1 by the nth expert is X ðnÞ q 1 �����q pÀ 1 ; Step 3: Constructing expert consistency networks for each factor The expert consistency network G Þ is formed when the group experts judge the influence degree of factor i on factor j in the subsystem F q 1 �����q pÀ 1 .The network adjacency matrix is calculated as Eqs ( 9) to (10); Step 4: Calculation of expert weight matrix by weighted network clustering coefficients According to Eqs ( 13) to ( 16), the clustering coefficient vector C q 1 �����q pÀ 1 ij of all experts for the network G q 1 �����q pÀ 1 ij ¼ ðV Þ is calculated, and the weights of all experts in this network are obtained by normalizing Eq (18); the consistency network for all factors is traversed, and the weight matrix W q 1 �����q pÀ 1 n of the experts is solved; Step 5: Calculate the weighted IDR matrix Calculate the weighted IDR matrix � X q 1 �����q pÀ 1 for subsystem F q 1 �����q pÀ 1 from Eq (20); Step 6: Solving the super IDR matrix by hierarchical DEMATEL method Using Eqs (5) to ( 6), the super IDR matrix is solved through weighted IDR matrix � X q 1 �����q pÀ 1 ; Step 7: Using the traditional DEMATEL method to calculate the centrality and causality of each factor According to Eqs (1) to (4), the centrality r i +d i of factor f i is calculated, which characterizes the relative importance of factor f i in the system; the causality r i −d i of factor f i is calculated, and if r i −d i >0, then f i is the cause factor, and if r i −d i <0, then f i is the receive factor.

Case presentation
There are many influence factors for combat capability, and the relationships between these factors are complex and intertwined.The hierarchical characteristics of these factors are obvious, so combat capability is a typical complex system.
Identifying and analyzing the key factors that influence combat capability is crucial for improving it.Sixteen factors (f 1 ~f16 ) that influence combat capability can be identified and classified into five dimensions: communication F 1�1 ¼ ff i ji ¼ 1; 2; 3g, intelligence F 1�2 ¼ ff i ji ¼ 4; 5g, command F 1�3 ¼ ff i ji ¼ 6; 7; 8; 9g, logistics F 2�1 ¼ ff i ji ¼ 10; 11; 12g, and fire support F 2�2 ¼ ff i ji ¼ 13; 14; 15; 16g, as shown in Table 2. Based on their attributes, these five dimensions can be considered as belonging to the command and control and communication systems F 1 ¼ fF 1�q 2 jq 2 ¼ 1; 2; 3g, as well as the fire and logistics support systems Communication, intelligence, command, logistics, and fire support factors are all critical in combat, and there are complex interrelationships among the factors within each subsystem, as well as between different subsystems.
The proposed method used to identifying the key factors influencing combat capability.In this case, combat capability includes 16 factors, and using the traditional DEMATEL method would require experts to make 240 judgments, which is obviously a significant workload.However, using the hierarchical DEMATEL method only requires experts to make 55 judgments, which clearly demonstrates the advantages of this method.
Step 1 Hierarchical decomposition The above elements are divided by hierarchy to form the structure diagram shown below The combat capability system is decomposed into a two-level structure according to the hierarchy, with the first level containing two subsystems, command and control and communications F 1 , and firepower and logistical support F 2 .The second level contains 5 subsystems of communication F 1�1 , intelligence F 1�2 , command F 1�3 , logistics F 2�1 , and fire support F 2�2 .These five subsystems specifically contain these 16 specific factors.Step 2 Constructing the IDR matrix for each subsystem The IDR matrix needs to invite experts to judge the influence relationship between the factors contained in the system, using the 0~4 scale method, 20 military theory researchers, weapon equipment professionals, combat commanders were invited to judge the system F,F 1 , F 2 ,F 1�1 ,F 1�2 ,F 1�3 ,F 2�1 ,F 2�2 , each expert needs to make 55 decisions, and the matrices X 1 , X ðnÞ 2 , X ðnÞ 1�1 , X ðnÞ 1�2 , X ðnÞ 1�3 , X ðnÞ 2�1 , X ðnÞ 2�2 respectively represent the IDR matrix obtained from the judgment of the nth expert on the corresponding system.
The decision situation of expert 1 is shown in Fig 7 .Due to the large volume of data, please see Appendix A at https://osf.io/gxtj5for additional expert decision information.
Step 3 Building expert consistency networks Taking the construction of expert consistency network G 11 as an example, the naming rule of the expert consistency network indicates that the superscript of G 11 represents the system number, and the subscript represents the factor number.Therefore, the consistency network formed by all experts when judging the degree of influence of the first factor of system F (i.e., system F 1 ) on itself.
The IDR matrix X (n) of 20 experts is extracted, and the element x ðnÞ 11 n = 1,2. ..,20 in the first row and first column is calculated according to Formula (9) to obtain the consistency matrix Y 11 formed by the 20 experts' scoring on this factor.The specific values corresponding to each row and column are shown in Table 3, and the expert consistency matrix for other elements can be found in Appendix B at https://osf.io/gxtj5.
In order to see the sparsity between experts more intuitively, according to the consistency matrix Y 11 of the network G 11 , we use Gephi to draw the structure of the consistency network G 11 as shown in Fig 8, even the thickness of the edge represents the weight of the edge, from the figure can be roughly seen that the experts are not equally close to each other.
Step 4 Calculation of expert weight matrix by weighted network clustering coefficients Continue to use network G 11 as an example to illustrate.The clustering coefficients of the 20 experts for the network G 11 are calculated according to Eqs (13) to ( 16) in Table 4.
The weights of the experts are obtained by normalizing the clustering coefficients according to Eq (17), the normalization index m = 20, to serve the purpose of reducing the weights of the experts who are far from the group consensus and giving more weights to the experts with high consensus.As can be seen from Table 5, the 8th, 11th, 15th, and 20th experts have significantly lower weight values than the other experts, which means that when judging the degree of influence on the system F 1 itself, the opinions of these experts are clearly inconsistent with other experts, and this situation can also be found from the matrix in Table 3, where the expert 20 with the lowest weight value, for example, he has a consensus degree of 0.75 with only 3 experts and consensus degree with other experts degree are all only 0.5, resulting in his low weight when scoring this element.
Due to the limitation of space, detailed data for the remaining expert weighting matrices can be found in Appendix C at https://osf.io/gxtj5.
Step 5 Calculate the weighted IDR matrix According to Eq (20), the weighted IDR matrix of each system is obtained as shown in Fig 10.
Step 6 Solving the super IDR matrix using hierarchical DEMATEL method Combined with the weighted IDR matrix of each subsystem, the super IDR matrix can be calculated according to the hierarchical DEMATEL method and Eqs (5) to (6) in Section Hierarchical DEMATEL method.The super IDR matrix integrates the mutual influence relationships of all factors, and the calculation results are retained to three decimal places, as shown in Table 6.
Step 7 Bringing the super IDR matrix into the traditional DEMATEL method, using Eqs (1) to (4) can be calculated to obtain the reasonability and centrality of each factor, and non-Archimedean infinitesimal is ε = 0.00001.As can be seen from Table 7, factors f 4 , f 6 , f 9 , f 13 , f 14 , f 15 , and f 16 are receive factors, and the remaining factors are cause factors, which is consistent with our common sense.
Table 8 shows the centrality of each element, based on which we can further obtain the weights of each element.According to the weights in Table 9, it can be concluded that the importance ranking of factors a influencing combat capability is It can be seen that the commander's command ability, communication interference and anti-jamming ability, the degree of intelligence of the accusation system, the adequacy of the war reserve, and the speed of attack are the most critical five factors.

Comparative analysis
Algorithm stability and deviation introduction.In most studies related to decisionmaking methods, scholars often analyze the superiority of the proposed method from a subjective perspective by explaining the reasons for differences in decision-making results.However, as decision-making problems themselves are subjective judgment problems, analyzing differences in decision-making results from a subjective perspective cannot provide sufficient evidence for judgment, as the analyst tends to analyze in their favor.This paper provides a new approach by comparing the robustness of decision-making methods in interference environments to analyze the advantages and disadvantages of the methods, i.e., comparing the stability and deviation of decision-making results in the interference of original data.The interference of original data is used to simulate the problem of random deviation in expert scoring caused by interference factors.Although this random deviation is a low-probability event, it can still cause changes in decision-making results when it occurs.Therefore, it is necessary to consider the ability of decision-making algorithms to maintain stability and the degree of deviation towards the true value.
The detailed definitions of the stability and deviation of decision-making algorithms are: The decision result of the complex system F ¼ ff i ji ¼ 1; 2; . . .; mg is obtained in descending order of weights ψ ¼ ½ f a 1 f a 2 . . .f a m �, and the subscript number of each factor corresponds to the decision vector η ¼ ½ a 1 a 2 . . .a m �.The decision result obtained from the original data in the interference-free condition is denoted as ψ 0 and the decision vector is η 0 .
Suppose that L decision experiments are conducted in the interference scenario, and the L decision experiments are independent of each other, and the decision results are ψ 1 ,ψ 2 ,. ..ψL , and the corresponding decision vectors are η 1 ,η 2 ,. ..,ηL .
Definition 1 Decision algorithm stability The cosine of the angle of the decision vector η 1 ,η 2 ,. ..,ηL in the Euclidean space is considered as the degree of agreement between the two decision outcomes in the L decision experiments, and for the decision vectors η i and η j , the degree of agreement between the two vectors is When cosθ ij = 1, it means that the decision vector η i is completely consistent with η j and the degree of consistency between the decision vectors has symmetry, i.e., cosθ ij = cosθ ji .From this, we can obtain the consistency matrix of L experiments The stability of the algorithm for L experiments is The matrix I is a matrix with all elements of 1, I2R L×L , denotes the most ideal case.When P = I, it means that all the decision results are equal when L experiments are conducted, and the stability of the algorithm is infinite at this time, which means that the decision algorithm is highly stable and the decision results do not change in any random interference scenario.
Definition 2 Deviation of decision-making algorithms The decision vector for the i th experiment in L decision experiments is η i , then we have b ¼ β is the deviation of decision-making algorithms in L decision experiments, which indicates the cumulative deviation of the experimental results from the true value.
As can be seen from the mathematical expressions of stability and deviation, stability mainly describes the degree to which the algorithm keeps the decision results consistent and stable with each other in multiple experiments.When the stability of the algorithm is high, it means that the decision results do not vary much from each other, indicating that the algorithm does not change the decision results easily due to random disturbances.The greater the deviation, the greater the deviation of the decision result from the true value after interference.
For the stability and deviation of a certain algorithm, there are several possibilities: 1. Large stability and large deviation Implying that although the results of multiple experiments are stable, they deviate far from the true value.

Large stability and small deviation
Meaning that the results of multiple experiments are stable and each experiment is close to the true value, which is the most desirable situation.

Small stability and large deviation
Meaning that the results of multiple experiments are not stable, and the results deviate from the true value, which is the worst case.4. Small stability and small deviation Meaning that although the combined results are close to the true value, the algorithm is not stable enough and is susceptible to random interference.

Comparative analysis of different methods.
To analyze the scientific validity of the method proposed in this paper and its superiority compared to other methods, we compare the fixed-point iteration method for calculating expert weights in reference [36] with the method proposed in this paper.This part continues to use the data in Section Case presentation as the original data, adding data perturbation scenarios for multiple experiments, each experiment is independent of each other, and the data used in both methods in the experiments are guaranteed to be exactly the same after the perturbation.
When the perturbation is zero, the decision results calculated according to the fixedmethod of immobile points are Some differences can be seen with the decision results calculated by the method in this paper, however, we cannot compare the advantages and disadvantages of these two methods by subjective analysis of the decision results only.Therefore, in this section, the analysis is based on the stability and deviation.
We set the perturbation environment as follows: The extent to which the data are perturbed is controlled by two parameters p 1 and p 2 , p 1 2[0,1] denotes the proportion of experts that are perturbed and p 2 2[0,1] denotes the proportion of elements that are perturbed in perturbed experts.As an example, we set p 1 = 0.3, p 2 = 0.2, p 1 means that 30% of the experts will be disturbed, p 2 means 20% of these experts' decision data will be randomly perturbed.The perturbed data will be randomly updated to an integer in the range [0-4], and the probability distribution follows a uniform distribution.
In order to consider all cases as much as possible, a range of p 1 20.04~0.36,p 2 20.04~0.36 with 0.02 steps was set and a total of 281 sets of experiments were performed considering different combinations of the two parameters.40 randomized experiments were conducted in each group, and the averaging method, the fixed-point iteration method and the method of this paper were compared and analyzed with the same data and conditions.The averaging method is the control group, and this method treats all experts' weights as equal, and in this paper, each expert's weight is set to 0.05.
From Table 10, it can be seen that inf appears in the second row, indicating that the stability of the algorithm reaches infinity at this point.According to Formula (22), when the results of L independent experiments are exactly the same, matrix P = I, which means that the stability of the algorithm is infinite, as indicated by Formula (23), and the corresponding offset of the algorithm is 0. This means that throughout the L experimental process, all experiments are exactly the same as the original value.
The first 10 sets of experimental data are given in Table 10, where method 1 refers to the averaging method, method 2 refers to the method in this paper, and method 3 refers to the method in [36].Due to the large amount of data, the detailed data of 281 sets of experiments regarding stability and deviation are recorded in Appendix D at https://osf.io/gxtj5.According to the results of 281 sets of experiments, the methods with the largest stability and the smallest deviation in each group of experiments were recorded separately, and the number and percentage of occurrences of the three methods were counted, as shown in Table 11.
As can be seen from Table 11, in the comparison of 281 sets of interference experiments, the method proposed in this paper has the largest stability value 199 times, occupying 68.86% of the number of experiments, and the smallest deviation value 236 times, accounting for 81.66% of the total number of experiments.It means that among these three methods, the method of this paper has the highest stability and accuracy, the method in [36] is the second, and the method of averaging directly on the expert decision matrix has the worst stability and accuracy.

Conclusion
In this study, we proposed a group hierarchical DEMATLE method for the identification of key factors of complex systems.The method inherits the advantages of the hierarchical DEMATLEL method, which can effectively reduce the workload of experts and, at the same time, large-scale group decision making enables more scientific and comprehensive decision results, which involves the expert weight matrix solving method to bring new ideas for the weight calculation when group experts make decisions.The main contributions and innovations of this paper are as follows: 1. Taking into account the expertise and limited knowledge of experts, the experts are assigned weights by factors to measure the overall performance of experts more finely with the weight matrix.2. The consensus of group experts is described by constructing an expert consistency network, and the degree of consensus of experts is expressed by the assigned clustering coefficients as an important basis for calculating the weights.
3. Stability and deviation indexes are proposed to test the effectiveness of the algorithm, which makes the decision algorithm test more convincing.
In this paper, the proposed method is applied to the identification of key factors of combat capability complex systems, and the proposed method is compared with other methods, and the experimental results have achieved good results.However, this article does not take into account the adjustment of expert opinions, and lacks the decision adjustment and opinion correction process of experts in reaching group consensus.This will be the focus of future research.